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SUMMARY 


The report covers the results of an analytical and experimental investi- 
gation on the vibrational energy transfer between connected substructures 
under random excitation. In the analytical area, the basic foundation and 
assumptions of the Statistical Energy Analysis (SEA) method, a major tool 
in random response analysis of structures, were examined and reviewed. A 
new SEA formulation based on the strong coupling condition of the substruc- 
tures was carried out and presented. Also presented were the results of 
vibration energy transfer study based on the wave equations applied to 
connected structures. In the experimental phase, three simple structural 
models were fabricated and tested. Additional tests were performed on 
selected substructures which formed parts of the test models. The test 
results were presented and evaluated against the analytical data. 

The work described in this report was carried out at Northrop Corporation 
under the sponsorship of NASA Marshall Space Flight Center. The contract 
number was NAS8-28171. The program was monitored by Drs. Hugo Steiner 
and Rudolph Glaser of NASA, under the overall direction of Mr. Richard 
Schock. The experimental work was carried out by Mr. Paul Finwall of 
Northrop. Dr. Paul Seide participated in the analytical work which is 
described in the Appendices of the report. 
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INTRODUCTION 


This report covers the analytical and experimental results obtained 
in the conduction of the Contract NAS8-28171 entitled "Investigation of 
Vibrational Energy Transfer in Connected Structures." The purposes of the 
subject program are (1) to investigate the high-frequency energy transfer 
mechanism between two connected structures, (2) to determine the major 
parameters that affect the energy transfer, and (3) to determine the para- 
meters and their influences to the application of the statistical energy 
analysis (SEA) method to the structures. To accomplish the program tasks, 
coordinated experimental and analytical methods were applied. Specifically, 
conditions were established under which the SEA method can be used to predict 
the vibrational energy transmission in two connected structures. Guidelines 
were established. A new analytical formulation was introduced to cover the 
cases where the substructure intercoupling was strong. In addition, analyti- 
cal and experimental definitions of the key parameters relating to the appli- 
cation of SEA were examined and applied to structures under high-intensity 
acoustic or random mechanical loads. 

The statistical energy analysis (SEA) method was first applied to the 
structural vibration problem as an extension of the room acoustics approach 
in acoustic engineering. Developed in the past decade by Lyon, Smith, Dyer 
and associates (References 1-15), the method considers the linear responses 
of multimodal structures and the resulting energy flow among the modes of 
two or more sets of substructures. The modes of a substructure are called 
a subset. Because of the complexity of the high-frequency modal data of 
structures and the desire to have a simple tool for engineering application, 
a number of assumptions are made in the formulation of the SEA. These 
assumptions, together with the related parameters which affect the validity 
of the assumptions, are considered to be the key points in determining 
the applicability of the SEA method to high-frequency vibrational energy 
transfer problems for connected structures. 
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In addition to the fundamental derivations of Lyon and associates 
described above, Zeman and Bogdanoff presented an elaborate formulation of 
SEA using technical terms and analytical procedures commonly used by struc- 
tural dynamicists (Reference 16). Zeman 1 s derivation deviates from the 
original SEA formulation in some details relating to the coefficients of 
dissipation functions, etc. While the end results are essentially the same 
as compared to the Lyon and Maidanik derivation (References 1, 11), Zeman 1 s 
derivation and its reasoning are helpful in making the SEA method more com- 
prehensible to the structural dynamicists. 

The SEA is based on the power flow between groups of linear oscillators 
Between two groups, the power flow is established by a set of dynamic 
equations. Each equation represents the mode response of one oscillator 
and its weak coupling with one or more oscillators from the other group. 

The coupling parameters are classified into inertia, damping, and spring 
types. For a stationary process, the assumption that the damping coupling 
parameters for any two oscillators are equal in magnitude and opposite in 
sign gives rise to a condition called gyroscopic coupling. Specifically, 
a gyroscopic coupling element is defined as one which produces a negative 
coupling force on Oscillator No. 2 due to a positive velocity of Oscillator 
No. 1 if the velocity results in a positive force on Oscillator No. 1 due to 
a positive velocity of Oscillator No. 2 (Reference 11). The gyroscopic 
coupling force on Oscillator No. 2 due to a positive velocity of Oscil- 
lator No. 1 if it results in a positive force on Oscillator No. 1 due to 

of the power flow coefficient under the weak gyroscopic coupling condition. 

Consider a narrow frequency band for which the modal density of the 
substructure may be determined either experimentally or analytically. In 
SEA, it is assumed that the input power spectrum is fairly flat within the 
frequency band. Each linear oscillator which is directly excited by the 
external source is considered to be subject to a "thermal bath." Under 
this condition, the modal energies of all the oscillators whose natural 
frequencies lie within the narrow band are fairly equal and may be repre- 
sented by an average value. A final formulation of the SEA involves the 
response levels of two or more substructures (which may be either connected 
substructures or a structure and a reverberant acoustic field) based on the 
average modal energies of the externally excited and the coupled oscillators 
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The SEA method has the potential of being a powerful engineering tool 
because of its generality and, specifically, its averaging technique, where 
the structural and response details are considered only in a broad sense. 

The method may also be misused if applied indiscriminately. In view of this 
background, the subject program investigates the key factors affecting the 
SEA, and specially the high frequency vibration energy transmission, so that 
the applicability of the method may be clearly defined. Furthermore, new 
analytical derivations are established which tend to clarify certain aspects 
relating to the formulation and application of the SEA. 

As part of the engineering method development, a user-oriented 
preliminary test procedure is developed in the application of SEA. The 
purpose of the preliminary procedure is to ensure that the structural model 
and the substructure definition implemented by the user will satisfy the 
basic assumptions of SEA. The same procedure also yields guideline indica- 
tions when the limits of the application of SEA have been surpassed due to 
such factors as substructure design, the operating frequency range, and the 
like. 

In the experimental phase of the investigation, three structural models 
were fabricated and tested in order to extract the maximum amount of in- 
formation from the test program. The models were designed to consist of 
connected substructures with typical variations in interface configurations, 
modal density distributions, vibrational energy transmission paths, and 
degrees of modal energy diffusion. 

The remaining text of this report follows the logical sequence of 
development. After a general discussion, the basic formulation of SEA is 
presented where the degree of coupling of the connected structural sets is 
discussed. This is followed by a new derivation on the SEA applicable to 
the strong coupling case. The subsequent sections cover the user guide- 
lines for SEA and the analytical formulations involving various types of 
connected structures. The text is concluded with the presentation and 
discussion of the experimental data. 
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DISCUSSION OF THE PROBLEM 


The current investigation emphasizes a detailed understanding of the 
mechanism of energy transmission at the interface of connected structures. 

In general, the high-frequency energy transmission is in the form of bending 
and shear waves for thin-skinned space vehicle structures. The tension/ 
compression waves may be involved in the picture, but they are considered 
to contribute minimally to energy transfer. The major parameters affecting 
the magnitudes and modes of energy transfer include the makeup of the sub- 
structures, their boundary conditions, the interface configurations (length, 
geometry, method of fabrication, etc.), the substructure and coupling loss 
factors, the relative amplitudes of the modal densities of the connected 
substructures, and the location and type of loading. In order to sort out 
the various parameters and to reach a rational solution of the complex 
problem, a number of simple test models featuring certain basic similarities 
were used. The details of the models will be described in a later section. 

In a given connected structure, the degree of modal diffusion depends 
greatly on the wave length relative to the characteristics dimensions of 
the substructure and, for a thin-skinned substructure, on its thickness. 

As the stress waves propagate in the structure, they meet the structural 
boundaries and interfaces where the waves are partially or totally reflected 
according to the boundary geometry and the constraint conditions. The 
infinitely many possibilities of these wave propagations and reflections 
cause a randomly distributed wave pattern (i.e., a high degree of modal 
diffusion). For waves of medium or long length, the degree of diffusion 
will affect the energy transfer through an interface because the directional 
properties of the waves determine the amount of energy transmitted to the 
neighboring structure. The effect is believed to be less pronounced for 
shorter- length waves. 

To confirm the energy transfer mechanism in a structural interface, 
analytical methods are applied using the classical equations of wave 
propagation. In a previous work by Lyon and Eichler (Reference 4), a set 
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of flexural wave equations was established. Simplified boundary conditions 
were applied to the interface. The energy transfer rate through the inter- 
face was estimated based on the assumption of thoroughly diffused waves. 

In the current program, more elaborate equations (e.g., those given by 
Mindlin in Reference 50) and more complicated interface configurations are 
used. For shell-type structures, the curvature effect and the proper con- 
straint at the interface are considered. Typical applications of the wave 
equations are described later in the report. 

In previous work on the application of SEA, it was found that the weak 
coupling conditions are not always satisfied, depending on the frequency 
range and other pertinent factors (Reference 49). In addition to estab- 
lishing the applicability of the SEA method for specific structures, alter- 
native approaches are investigated and reported where non-weak intercoup- 
lings of subs true tures are involved. 

To make the SEA method a usable tool for practicing engineers dealing 
with high-frequency vibrations of connected aerospace structures, it is 
desirable to have available general guidelines on the proper application of 
the method. The guidelines presented in this report take into consideration 
such structural parameters as coupling loss factors, damping loss factors, 
modal densities, etc. These guidelines also include preliminary test pro- 
cedures which can be followed by practicing engineers. In the subsequent 

sections, the basic foundation of the SEA method, the applicability of the 
method to strongly coupled structures, and the guidelines and preliminary 
test procedures to the proper application of the SEA method are described 
in detail. 
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BASIC ASSUMPTIONS OF THE SEA METHOD 


The theoretical foundation of the statistical energy approach to vibra- 
tion analysis is that the steady-state time-average power flow from one mode 
to another is proportional to the difference between the time-average kinetic 
energies of the two modes provided that the following conditions are satisfied 
(References 1, 11): 
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The coupling between the two modes which satisfy the following set 
of equations is linear, weak, and conservative (gyroscopic). 


ra l (y l + \ + w l ' 2 y l> + A y 2 + B 2 y 2 + c y 2 = f l 

m 2 ^ y 2 + ^ 2 y 2 + *2 y 2^ + A y i + B i y i + C y l = f 2 


( 1 ) 


where m, ( 3 9(j j are the mas s, damping coefficient, and the natural fre- 
quency of the oscillators, respectively. The modal displacements are 
denoted by y. In Equation (1), the coupling parameters A and C rep- 
resent the inertial and stiffness couplings, respectively. The coupl- 
ing parameters B^ and represent the gyroscopic coupling when 

B i = "V 

2. The forces f^ and f 2 (see Equation (1)) acting on the two modes are 
uncorrelated and have spectra that are relatively flat within the 
frequency band encompassed by the resonances of the coupled system. 


Furthermore, the steady-state time-average power flow from one set of 
modes to another is proportional to the difference between the set average 
modal kinetic energies of the two sets provided that either: 

1. The mode-to-mode coupling is the same for all mode pairs, or, 

2. All modes in each set have equal time-average kinetic energies. 
(Note that the members of each set are not coupled to each other.) 
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The coupling factor, applicable to set-to-set power flow is equal to the 
sum of all the mode-to-mode coupling factors. 

For a given connected structure, the SEA method can be applied to predict 
the response levels under high frequency excitation. For this case, it is 
further required that the wave patterns are diffused, and the major wave 
lengths of interest are small compared with the characteristic dimension 
of the structure. Based on the above, the prerequisite conditions of appli- 
cability of the SEA to a connected structure may be restated in terms of the 
analytical formulations given below: 


1. The generalized coordinates of displacements of the connected 
structural system satisfy the following set of equations: 

( -jUj + Pjij + ^ 2 x t ) +J [a,^ + + C k .y k ] = fj 1-1, 


, . . . * , N 


'j (y j + + "j 2y ) ) + B /1 i l + f i 3 = 1 8 


where A., = A, , and C,, = C , . 

Ji ij ji 1 J 

The coupling is gyroscopic, i.e., B + B = 0. (3) 

The coupling is weak, | A^_. | , | A ^ ^ | « m^, m^ ; 

I B tj | , I B j± I «m 1 P 1 , m j P J .jC i J , \c ^ \«m i ^ i 2 , m jWj 2 . (4) 


In addition to conditions noted in Equations (2), (3), and (4), other 
conditions such as uncorrelated modal forces, equal modal energies, etc., 
are also needed in order to apply the SEA to the connected structure. 

In general, for a given connected structure, it is not always clear 
whether the values of the coupling parameters A, B, and C satisfy condi- 
tions contained in Equations (2), (3), and (4). 

In the following, the degree of coupling of a connected system used by 
Lyon and Eichler is examined. 
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In Reference 4, Lyon and Eichler treated analytically the coupled system 
shown in Figure 1, namely, a mass-spring resonator attached to a finite thin 
plate. The same configuration is used here to examine the magnitude of the 
coupling parameters of the coupled system. The transverse displacement w 
of the plate satisfies the following equation: 

Kp a Cp a V* w + w + ftp w = -i [r(x, t) - f(t) 6U - Xq)] 

r P s ( 5) 


where K p is the radius of gyration of the cross section, C p is the longi- 
tudinal wave velocity in the thin plate, P p is the damping coefficient 
related to the loss factor r| p by (3 p = n p P g is the mass per unit area, 
r(x, t) represents the random loads per unit area, and f is the reaction 
force produced by the resonator which is attached to the plate. The cor- 
responding modal response equations are: 

+ ^>m + ^m + ^ = 

Ps 

where = modal amplitude of m“ t ^ 1 mode 

- th 

U) m = natural frequency of m mode 


Y m ° = mode shape at x = Xq 

F m = modal force corresponding to m" ^ mode 


Using the reaction force as the dependent variable, the equation of 

motion for the resonator of Figure 1 is: 

• • • 

a f - K L Y m ° UB 3 f s (t) 


f + 9o f 


Wo 


m 


(7) 


r 2 K 

where 6 = — and to = — . The exciting force on the resonator mass is 

r o M o M 

f g (t). The reaction force f is related to the displacement of resonator 

d by the following relation: 

K(d - w 0 ) = -f = K(d - £ Y 0 g ) 

m n> m 


( 8 ) 


The corresponding power balance equations for the system are: 

< Ps^m >+ ^m ^ ^m > = Ps < ^m §m > 

< f a /K > -E V < f g > = I < fjfdt > *, - < f s d > 

* m m M 00 & 


( 9 ) 
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The above equations show that the coupling of the system is conservative. 
Since the equation of motion for the resonator is in a special form, Lyon and 
Eichler introduced the following new variables: 


e = 

(K 

Ps)" k 

f 

v m = 

(K 

Ps>* 


= 

(K/ 

P.> % 



Based on the above, a new set of equations are obtained as follows: 

v m + ftp v m * ® = ^ Ps^ F m ( t) (11) 

9 + fc 6 + ub 3 8 - L v m = u^ 2 6 s (t) 

Comparing with Equation (1), Equation (11) satisfies the following conditions: 
^im = ^mi = 0 

Bxm = -B mx = * m ° (12 ) 

m = Cmi = 0 


The above conditions establish the gyroscopic coupling of the system 
and the applicability of the SEA method. Based on Equation (11) and the form- 
ulas given in Reference 1, the amount of energy transfer through the spring 
resonator and the finite thin plate may be computed. Consider a simply sup- 
ported square plate with a resonator attached at the center. We have: 


*m° 


s ^ . B x m - 




(13) 


where a is the edge length of the plate and Mp is the mass of the plate. 

In the narrow frequency band with center frequency the ratio of 

B/pp is: 

B_ _ 2_ j M_ (14) 

Pp " \ V Mp 
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Figure 9 of Reference 4 shows r\ - 0.0005 for a 0.145 cm (0.057 n ) steel plate. 
So we have: 


B_ 


4000 



(15) 


It is noted that for the case described above, ra^ of Condition (4) is 
equal to unity. In other words, the weak coupling Condition (4) is now 
B/Pp << 1, which may be restated as: M«Mp/(16 x 10^). Since the latter 
condition is not satisfied in a typical set-up, we conclude that the resonator- 
plate arrangement is not a weak coupling case. In the following section, new 
formulation of vibrational energy transfer involving strong coupling of sub- 
structures is described. 



FIGURE 1 DIAGRAM OF PLATE WITH ATTACHED RESONATOR 
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THE EXTENSION OF THE SEA METHOD TO STRONGLY COUPLED STRUCTURES 


A major item of interest in the application of the Statistical Energy 
Analysis is the degree of coupling of connected structures* Previous work 
at Northrop indicated that the weak coupling condition assumed by the original 
workers of SEA was not always satisfied in typical structures. Our present 
effort is to determine the extent of coupling of practical aerospace struc- 
tures and to explore the additional formulation in structural responses based 
on the condition that the substructure coupling is not weak. Our findings 
are described below. 


Power Flow between Two Modes 

The average net power transferred between two linearly strongly coupled 
oscillators which satisfy the set of Equations (1) is derived below. The mean 
value function equations of Reference 1 are the starting point of this deri- 
vation. Thus, the sources f^ and f 2 are assumed to be statistically inde- 
pendent and to have power spectra which are flat as compared to the admittance 
spectra of the two oscillators. Suppose an impulse n a M occurs from source f^ 
only, then by Equation (1) 

f m.Ay. + A^y = a, 

11 (16) 

( AAy^ + m 2 4y 2 = 0. 


Accordingly, 



am 2 /(m 1 m 2 - A ) 
2 

-aA / (m^m 2 - A ) 


(17) 


From the nature of f^, subsequent increments of Ay ^ and^y 2 (due to 
different impulses) are independent increments of velocity, and therefore 
in one second the amounts of energy gained by oscillator 1, 2 are 

2 I 2 2 2 

m l < ^1 | = m l m 2 D i/( m i m 2 " A ) 

1 sec. (18) 

2 I o 9 2 

m 2 < ^2 | = m 2 A Dj/(n»jm2 - A ) 

1 sec. 

where is the spectral density of 
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Thus, the power from source 1 is 


<fjY^> = D 1 m 2 (m 1 m 2 4- A^/Cm^ - A 2 ) 2 (19) 

and similarly, 

<f 2 Y 2 > = D 2 m 1 (ra 1 m 2 + A 2 )/(mjm 2 - A 2 ) 2 (20) 

also, since 4y^ = one has 

<f 1 S r 2 > = -D 1 A(m 1 m 2 + A 2 )/^^ - A 2 ) 2 (21) 

and 

<f 2 yj> = -D 2 A(m 1 m 2 + A 2 )/(mjm 2 - A 2 ) 2 (22) 

Since an impulse produces no immediate increment of displacement, 

<f 1 y 1 > = <f 2 y 2 > ~ <f l y 2 > = <f 2 y l > = 0 

The above covariances between the sources, displacements, and velocities 
define the statistical properties of the sources. The effective force which 
the motion of oscillator 1 causes to be produced on oscillator 2 is seen from 
equations (1) 


f 12 = -^1 - Vl - Cy l 


(23) 


Thus, for stationary ergodic process the average net power flow delivered from 
oscillator 1 to mode 2 is 


P 12 = <f 12 y 2 > = - A< y 1 y 2 > ~ B 1 < y l y 2 > “ C <y l y 2 > 


( 24 ) 


Similarly, the net power transferred from oscillator 2 to 1 


P 21 =<f 2lV = ^ <y l y 2 > _ B 2 <y l y 2 >_ C <y 2 y l > 


For a stationary process 


and therefore we may write 

P 2i = A< y’iy 2 > ~ b 2 <y iV + c <y i^2 > < 26 ) 


In order to evaluate the moments involved in Equations (24) and (26) the follow- 
ing set of equations may be obtained from stochastic equations (1) 



Equation (27) may be solved explicitly and the results used to obtain power 
flow etc# This is accomplished if the coupling is conservative (gyro- 
scopic coupling = -B^ = B) . 

P 12 = ~ P 21 = g 12 (e l “ e 2 ) (28) 


where P 


ij 


G. 

l 


the average net power delivered from oscillator i to oscillator j. 
the stored energy in system i = <f^y^> / 


gl2 — g2i ~ the power flow coefficient 

= (m 1 m 2 P 1 2 p 2 2 | A 2 [Pl w 2 4 + P 2 W 1 4 + PiP 2 ( ^1 W 2 2 + ^ 2 2to l 2 ^ 

+ (B 2 - 2AC)(p lW 2 2 + P^j 2 ) 

+ c 2 (P l + p 2 )} + b 2 P 1 P 2 {a 2 [p ^ 4 + P ^ 4 - P 1 P 2 (P 1 co 2 2 

+ B^ 2 )] - (B 2 + 2ac)(b iW2 2 + e^ 2 ) + C 2 (gj + 6 2 )| ) . 
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i m l m 2 B l 6 2 (w l 2 " w 2 2)2 + (>» 1 m 2 6 1 B 2 - B 2 - 2AC) • 

<0 ^ + B 2 ) (Bj^ 2 + + + ^ 2 ^ C + 

(B 1 u 2 2 + B^ 2 ) 2 A 2 } -1 (B 2 4tn 1 m 2 B 1 6 2 ) -1 (29) 

For the special case of weak coupling between two oscillators, 

2 2 

( |A|«m 1 ,ra 2 ; |C|« m^w^ , ; Ib m 2^2 and m l = ^ n 2 S==1 ^ • Equation (29) 

may be reduced to the form which is identical to the result of Reference 1, 

Equation (3,4) • 

Equation (28) shows the fact that the steady-state time-average power flow 
from one mode to another is proportional to the difference between the time- 
average modal energies even though the coupling between the two modes is 
strong, provided that all other requirements to apply the SEA method are satisfied. 

In the following subsection, the power flow between two strongly coupled 
sets of modes is described. 

Power Flow between Two Strongly Coupled Sets of Modes 

Consider two sets of modes. Within each set, the modes are uncoupled 
to each other. The power flows between modes of the two sets are assumed to 
be proportional to the modal energy difference. The modal displacements of 
the oscillators, which are denoted by x^, y respectively, satisfy Equations 
(2) and (3). Furthermore, assuming a stationary process, the time average 
of a function is denoted by a pair of brackets < > around the function. 

The power balance equations may then be expressed as follows: 


j m i fi i 

(*?) 

= P. e. - 

i i 

N 

g ik^ d ik 

- V 

i = 1, • • • • , N 

(30) 

i«r 

i a 


Vj + 

N 

-V 

j = 1, • • • • , N 

(31) 
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or 


MK)- e i%-Wu-V 

]*,{*■)= + g lj < e ’u - 


(32X 


i = .1, ■•••, N;j = 1, • • • , N (33) 

where 0^ , 0^ are called the equivalent modal energies as defined below: 


and 


= 9 i - E *ik <9 'ik - "id > / "i 
k*j 

9 j'i = *J . g t) ' *j/> /? j 


«i = (fjX,) //», 

? i "<V J> /? i 


(34) 

(35) 

(36X 

(37) 


Thus, the expression of the mode-to-mode power flow and the relation 
between coupling parameters and g are as shown below: 

<p ij > = 8 U <e ij - 5ji> = - A ij < Vj >_ B ij < Vj > ' °ij < Vj > (38) 

where the power flow coefficient g^ is defined by Equation (29). Equations 
(34), (35) and (38) are used for defining the energy terms. Based on this 
formulation, the power flow diagram for a typical mode (x^) and the mode (y ) 
is shown in Figure (2a). The steady-state time-average power flow from one 


set of modes to another may be expressed as 
N U 

~~ g ij ^ 6 i j ” 

i=l j=l 


IJ 


(39) 


Taking the summation of the terms of Equation (33) with respect to index 
"i" from 1 to N yields 

.2 v r v*' , . r» (40) 


Nm j^(^ 2 )= + t ? Kj . yj,) 
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Insertion of Equation (31) into Equation (40) gives the following 


expression 

N 

2 ^ + (» - D < 

i=i 


• 2 > 

y j > 


( 41 ) 


— I 


In Equations (40) and (41) the unknown 0. may be represented in terms of 

* 2*2 ^ 

<y^>and<x^>. This is accomplished by solving the simultaneous equations 

(32) and (33). 


e ij = t m j O- - 8ij/ 6 i) <^j 2 > “ < *i 2> J / A 


ij 


(42) 


i 2 2 

0 ji = t m l (1 " < *i > - < yj >] /Ajj (43) 


where 


A 


ij 



+ •i" ,) e ij 


(44) 


Substituting Equations (42) and (43) into Equation (39), the following expres- 
sion is reached 


P 


IJ 


N N 


Z E 

i=l j=l 



(45) 


It may be concluded from Equation (45) that the set-to-set power flow 
is proportional to the difference between the set-average modal energies 
of the two sets provided that either or both of the following two conditions 
are true: (1) the mode— to-mode coupling is the same for all mode pairs; 

(2) all modes in a set have equal time-average energies * Furthermore, Equa- 
tion (45) may be rewritten as 



» \ 

(— 

2 - 

-2 \ 


p = ( N 

A-* 8 ij^ A ij/ 
j=l J 

<V > _ M 

<V > 1 

(46) 

ij \^> 
1=1 

\ 

N 

N / 

where the M, N, 

2 

<V > are the mass. 

the 

mode count, 

and the mean 

square velocity 


— — —2 

of Set I; and M, N,<V > are the corresponding functions for the Set J. The 
energy flow diagram of two coupled sets is shown in Figure (2b). 
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FIGURE 2 ENERGY FLOW IN TWO COUPLED MODE SETS 


17 



In Reference 11, a term called coupling loss factor of Set I is defined 
to describe the loss of mechanical energy from Set I, the energy being trans- 
ferred to set J through the interface: 


n c = ^IJ 


IJ 


U) M <V > 

o 


N N 


—2 

< V > = 0 


^ N (Z) Z 8 ij /A tj) 

° '1=1 j=l 1 


(47) 


where w is the center circular frequency of the narrow band excitation. The 
° -2 

condition <V > = 0 in Equation (47) is imposed in order to make this loss 
factor independent of the modal energies of the system receiving energy 
inflow. 


It is noted that the coupling loss factor n = n TT of set J is generally 

^ J I 

not the same as the loss factor 


n = n T T of Set I. Instead, we have 
c I J 


N N 


N n ^ 
c o 


* n c"o =EE ®ij % 
1=1 j -1 


(48) 


Ways to determine the coupling loss factors analytically and experimentally 
for various connected structures will be discussed later in the report* 


Estimation of the Response Ratio of Two Connected Structures When One Is 
Directly Excited 

In the SEA method, an important application is to estimate the response 
ratio when the primary structure of the connected structural system is randomly 
excited in a narrow frequency band. The response ratio may be estimated by 
inserting equation (43) into Equation (41) and setting (Set J is excited 

through attachment only) : 


N 


E m i <x j > 




i=l 


N 


8 j + H g ij /A ij 

1=1 


(49) 
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The above result may be expressed in the following simplified form 


— 2 2 MN ® c 

< V > / < V > = jf 1 — — 

MN B + B 


( 50 ) 


where g is the average damping coefficient of the coupled structure J and 
6 is the coupling damping factor defined below: 


N 

B = mean of g../A,. for all j's 

c i=l 1J lj 


(51) 


^ r ' c w o 


c 

n JI w o 


Equation (50) conforms to the general formulation of SEA (see Equation 101, 
Reference 11). Furthermore, a term called the apparent loss factor n 
is introduced in Reference 11 (Equation 103). The apparent loss factor 
is the value of the dissipation loss factor ascribed to set I (on the basis 
of measurements performed on Set I) if the observer is not aware that the set 
is coupled to Set J, i.e., 

a n 7 t ^ n c ^ 

n = n + 1 = n + (52) 

J n JT + ^ |N j n c + n 

where n and n denote the loss factors of Set I and J respectively. It is 
evident from equation (52) that the apparent loss factor n T is always greater 

J- J 

than the actual dissipative loss factor 

Based on the above equation, the mean square of response ratio of the 
two mode sets (Equation (52)) may be reduced to: 


< V > / <V > 


M 

M 



(53) 
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Conservation of energy requires the following conditions. 


n „o«< v > +P ,J 


= n TT U) M <v > 

IJ o 

:v 2 > 


IJ 


= n n 


/ M <V > _ M< V > \ 

N N ' 


_2 

n M <v > 


(54a) 

(54b) 

(55a) 

(55b) 


where is the time-average power supplied to Set I (Figure 2b). Equation (55a) 
is an alternate form of Equation (46). 


In Reference 4, Lyon and Zichler applied Equation (50) to two coupled 
plates (Figure 3a) using the apparent loss factor: 


a 



(56) 


Comparing the above with Equation (52), it may be concluded that Equation (56) 
is true only when the weak coupling condition is satisfied. The weak coupling 
condition is: 


n 


ji 


n « n or n « JLri 
c U N 


(57) 


In Equations (56) and (57), Set I denotes the vertical plate while the set 
J denotes the base plate of the two— plate system. 

In the test phase of the program the loss factors and the apparent loss 
factor were measured experimentally for various structural models. The acquired 
data and the conclusion drawn from these data will be given in a later section. 
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GUIDELINES AND PRELIMINARY TEST PROCEDURES 


TO THE PROPER APPLICATION OF SEA METHOD 

In order to make the SEA method a usable tool dealing with high- 
frequency vibrations of connected aerospsce structures, it is desirable 
to have general guidelines on the proper application of the method. The 
guidelines described here include a preliminary test procedure which can 
be followed conveniently by practicing engineers. The purpose of the 
preliminary procedure is to ensure that the structural model and the sub- 
structure definition implemented by the user satisfies the basic assumptions 
of SEA. The same procedure also yields guideline indications when the limits 
of the application of SEA have been surpassed. The procedure is formulated 
below. 

The theoretical foundation of SEA in vibration analysis demonstrates 
that under certain conditions, the average rate of flow of energy between 
two sets of modes (representing groups of modes of two coupled systems in a 
given frequency band) is proportional to the difference in the set-average 
modal energies (Equation (46)) 


„ - - / M<V > 

p u = N N g u ( “ — 


M<V 2 

N 


(58) 


where denotes the average value of the power flow coefficients from a 

mode of set ”1” to a mode of set "J." Equating Equations (55b) and (58) 
yields 


<v 2 > = M_N N S IJ ( 59 > 

<V 2 > MN nw+Ng TT 
O IJ 


The above equation is equivalent to Equation (50) where "J" represents the 
substructure excited through the connecting interface only. Furthermore, 
by means of Equation (48) we have 

«U - \ V" (60> 
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For a connected structural system, if all the basic assumptions of SEA 
are satisfied, it may be concluded from the above equations that 

*ij**ji >o <6i) 

The proposed test procedure will check whether g (Set "I" excited) 
is approximately equal to g J]; (Set "J" excited) for all frequency bands of 
interest. (The ways to determing g values will be described later in the 
discussion.) It is expected that the test results will fall into one of 
the two categories given below: 

1. The values of g differ substantially from the values of g It 
indicates that the SEA method cannot be applied to the connected 
structural model under test. 

2. The values of g are approximately equal to the values of g . It 

t J j I 

indicates that the coupling between the connected structures is 
conservative and the SEA method may be applied. 

For typical structural elements, mass M and modal density n are known 
quantities. The dissipation loss factor H may be determined based on the 
measurement of the decay time. Another test method called the Q-method may 
be applied to determine "H, The total number of modes N in a set may be de- 
termined as the product of modal density n and the frequency band of exci- 
tation Ato. For a connected structural system after the substructure data 
M, N, T,M, N, *1 are determined, the values of g may be determined by either 
one of the following approaches: 

1. Excite primary structure "I" of the system at points chosen at 

2 —2 

random and measure the response levels < V >,< V >. The values 
of g will be computed based on Equation (59). 

2. Compute g^ based on equation (60). The values of the coupling loss 
factor r| c may be obtained as a function of *1, T, ^ i j ^y solving 
equation (52). 

Mi!. - n) 

\ " Z “J = (62) 

^(nij-WN 

The loss factor may be determined either experimentally or analyti- 
cally (Reference 4). 
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The aforementioned test procedures were performed on the three fabri- 
cated test specimens which are described in a later section of the report. 

The same procedures may be used to various types of connected structures* It 
will yield guideline information when the limits of the application of SEA 
have been surpassed due to such factors as structural element designs, the 
operating frequency range, etc. 

In the following section, analytical methods are applied to investigate 
the energy transfer mechanism in the structural interface of two types of 
connected structural systems. 
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VIBRATIONAL ENERGY TRANSFER MECHANISM IN 


TWO CONNECTED STRUCTURES 

In this section, the mechanism of energy transmission at the interface 
of two connected structures is investigated analytically. For this purpose, 
the classical equations of wave propagation are employed based on the assump- 
tion of thoroughly diffused waves. The results developed are compared with 
the test data given in a later section. Two typical structural models are 
considered : 

1. A two-plate system considering rotatory inertia and transverse shear. 

2. An infinite flat plate with an integrally attached half circular 
cylindrical shell. 

Two— Plate System Considering Rotatory Inertia and Transverse Shear 

The classical two dimensional theory of flexural motions of elastic plates 
is good only for waves which are long in comparison with the thickness of the 
plate. In case of transient loads with a sharp front, the significant fre- 
quencies of modes of vibration are of a high order. The flexural wave lengths 
of interest may reach the order of the plate thickness. An improved theory, 
which takes into account the effects of transverse shear deformation and 
rotatory inertia, should be used. In the following, the more elaborate plate 
equations of Mindlin (Reference 50) are used to determine analytically the 
energy transfer mechanism at the interface of a two-plate structure (Figure 
1-1) which was investigated by Lyon and Eichler (Reference 4) using the 
Bernoul 1 i— Eu ler plate equation. In the approach, the base plate is assumed 
to be continuous and supported along the interface in such a way that the nor- 
mal deflection vanishes and the twisting moment is continuous. The vertical 
plate hypothetically extends to the middle surface of the base plate where 
it is fastened along the interface. The detailed formulation of the flexural 
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wave propagation is given in Appendix I. If the effects of the transverse 
shear deformation and rotatory inertia are omitted, the result is identical 
to that obtained by Lyon and Eichler in Reference 4. 

In Reference 27, Heckl used an analogy of architectural acoustics to 
characterize the localized boundaries of plates in terms of absorption coeffi- 
cients. The absorption coefficient Y is defined as the ratio of the outgoing 
wave energy rate per unit length of the junction line vs. the incoming wave 
energy rate in the other plate. The absorption coefficient Y may be related 
to the difference between the apparent loss factor n and the loss factor n 
as shown below: 

YC L 
g 

(n a -n) w = (63) 

7T s 

where S is the area of the plate, L is the total length of the boundary, and 

C is the group velocity for bending waves on the plate, 
g 

In Appendix II, the effect of the coupling of the flexural and tangential 
waves on energy transmission is investigated analytically. In Appendix HI, 
the transmission of a flexural wave in an infinite plate with an integrally 
attached half cylindrical shell is analyzed. The end results of Appendices 
I through III are applied in Appendix IV where the energy transfers in various 
structural systems are formulated. Specifically, Appendix IV of the report 
shows the energy transfer mechanism in the two structural systems as previously 
mentioned. The expressions of the average input power and the transmitted 
power are formulated. It also shows generally the orthogonality relationships 
for the average power expressions. In other words, the average power involved 
in the forces of one mode of motion moving through the displacements of another 
mode is zero. 

For the two-plate system to be considered, as shown in Reference 4, the 
junction absorption coefficient may be computed as the total absorbed power 
averaged over all angles of incidence <f>^ (Figure I— 1)> and divided by the aver- 
aged incident power. In the following, the expression of Y for the incoming 
wave in the vertical plate is formulated and simplified. All the symbols used 
are defined in Appendix I. 
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Based on Equation (IV— 7) of Appendix IV, the average input power per 
unit width at an angle of incidence 4>- is given by 

r- * 

4 


P 1 = 


2 c 


— COS <t> 


1 


2 - \ 2 
C P 

1 1 


’1 


C l P l h l 


(64) 


The corresponding average transmitted power per unit width is given by 
Equation (IV-32) as 


P Q = 2 | | ^ cos <J >2 


4 

U) D, 


2 c. 


2 \ 2 4 

|l - li -Z | + C 2 P 2 h 2 


G 2 


a) 2 


(65) 


If sin 4> ^ = s * n ^1 > outgoing power vanishes. 

By definition, the absorption coefficient Y may be expressed in the following 
form 

J- tt/2 P o d< h 

(66a) 


Y = 


/ 


11 1 2 


-it/2 


p i d *i 


where 


H 2 P 2 h 2 C 2 Z * 71 / 2 
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h l P l h l C l ” 7T / 2 

1 - g t V 1 + 8 t + g t 2 


. + I 2 

A 2 I cos p 2 dij>^ 


(66b) 


n i - 

[Vl + g t 2 + s ( ] % 

(67) 

g i = 

I , 2 ) 

(68a) 

s i = 

1 *. (1 + L- v i * 2 ) 

2 1 2 

(68b) 

•i = 

P c* 2 / G ' 
i i ' i 

(68c) 
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where 



In general, for a thin plate with thickness h = 0(0*2 cm,) and the 

2 2 2 

interested frequency range f < 20 K Hz* The values of $ , g , s as de- 
fined in Equation (68) are much smaller than unity. Under this condition. 
Equations (69) and (70) may be simplified as follows: 

Y - 2r 


(1 - 2 gj) (1 - g 2 ) r V/2 cos'f'j cos <f> 2 d ^ 

(l + s 1 )(l+s 2 ) (l-g 1 )*' 7r/2 I Al 
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A 2 1 = 


a/1 + gi) 


(1 - Sj) 


^ a 2 d + 8 2 ) 

+ r 


n 2 


d - s 2 ) 


+ 


d - gj) 


v 


cos <f> 


1 


1 + s. 


1 - 


+ r 


yjl + s. 


cos <p 


(73) 


The expressions for a and may be represented as 


sin <t> 2 = * 


Vf 


+ S, 


sin 


+ s. 


(74) 


a 1 « V 7 


+ sin <t>- 


1 - 


s i cos <t> i 


2(1 + sin <(>j) 


(75) 


Two simple cases of Y of special interest to the present program are 
described below. The plates are assumed of the same material. 


1. h t = h 2 


2. h t « h 2 


(1 - 2g) ( 1 _ s) 

-1 + 

1 + a . - 

4> ( 1 + s ) 

sm 

vT 

h^/h^ <<: l f then 

3 (1 - 2bg )(1 - 

b 1 

- g 2 ) 

J 1 + s 2 

d-bg 2 ) ^l+s 2 +b s 2 

Y 1 + bs^ 


VT + 


(76) 
a = 2 $ /(l-s) 


(77) 


After the absorption coefficient is obtained, the corresponding loss factor 
data may be computed based on Equation (63). 
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Flat Plate with an Integrally Attached Half Circular Cylindrical Shell 


In the analytical work on the wave propagation between two perpendicular 
plates (Appendix I), the junction line has been assumed to remain straight 
and the coupling between flexural and extensional waves was ignored. Since 
flexural and extensional waves are always coupled in the cylindrical shell 
equations of motion, an analysis was performed where coupling between flex- 
ural and extensional waves is taken into account. The results are given in 
Appendix II. The formulation indicates that there is a range of angles of 
flexural wave incidence for which coupling effects are important. These results 
also indicate that the proper boundary conditions should be used for the plate- 
cylinder structural systems. 

The detailed formulation of the wave propagation in an infinite plate to 
which an infinite half circular cylindrical shell is integrally joined along 
the diametrically opposite generators (Figure III-l) is given in Appendix III. 

A harmonic flexural wave is assumed to be produced in the plate with the angle 
of incidence <J> which is scattered by the half cylindrical shell. For this 
case, in order to simplify the problem, rotatory inertia and transverse shear 
deformations are neglected in both the plate and the cylindrical shell. The 
proper boundary conditions, in which the coupling between flexural and exten- 
sional waves is taken into consideration, are used. 

The expression of the average input power per unit width and the average 
transmitted power in the half cylindrical shell are defined in Appendix IV, 
by Equations (IV-4C) and (IV-34) respectively. Some analytical computation 
based on the formulation of Appendices III and IV is given in the following 
section dealing with experimental results. 
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EXPERIMENTS ON VARIOUS STRUCTURAL SYSTEMS 


In the current program, a combined experimental and analytical investi- 
gation was performed on the random vibration and energy transfer in connec- 
ted structures* In the following, the experimental results are described. 

The description starts with the model fabrication and testing procedures. 

The test data acquired from the different models are illustrated and com- 
pared with the analytical results. 

Model Fabrication and Testing Procedures 

Three simply connected structural models and three separate components 
were fabricated in the performance of this experimental investigation. 

The specimens were fabricated from Gage 16 (.16 cm) and Gage 22 (.081 cm) Type 
304 stainless steel plate stock. Electron beam welding techniques were 
employed to fabricate the simply connected structures to minimize warping 
of the specimens. The first model has an irregular shaped vertical plate 
welded to a rectangular base plate. The second model has a half cylindri- 
cal shell welded lengthwise to the base plate along the longitudinal bound- 
aries of the shell. The third model consists of welding an open— end circu- 
lar cylindrical shell to the base plate. All three models were supported at 
the corners of the 16-gage rectangular base plate. The three components 

fabricated were duplicates of the model No. 1 vertical plate, Model No. 2 
half cylindrical shell, and a base plate which was common to all three 
models. The three models are illustrated in Figure 3 and photographs of 
the models are shown in Figures 4-6. Also shown in Figure 4 is a typical 
deformation pattern (center frequency 2668 Hz, half wave length 1.5 inches) 
yielded by the visualization technique through the use of lightweight poly- 
vinyl chloride particles. 

For each of the simply connected structural models and the components, 
the following test procedure was employed. A Goodman 390A shaker was attached 
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(a) MODEL NO. 1 

DIMENSIONS ARE TYPICAL FOR ALL THREE MODELS 


HALF CYLINDRICAL SHELL 
THICKNESS .079 CM (.031") 




FIGURE 3 SKETCHES OF TEST MODELS 
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FIGURE 



4 TEST MODEL NO. 1 AND THE TYPICAL DEFORMATION 

PATTERN OF THE BASE PLATE WITH CENTER FREQUENCY 
f = 2668 Hz, HALF WAVE LENGTH = 3.8 cm. 
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FIGURE 6 TEST MODEL NO. 3 
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through a Kisstler 931A force transducer to a selected point on the model. 
Endevco Model 2222B accelerometers were then cemented to the surfaces of 
the test specimen at selected locations. The specimen was then excited in 
1/3 octave bands ranging in center frequency from 250 Hz to 16K Hz. During 
random excitation in each 1/3 octave band, the rms response acceleration was 
recorded for selected locations on the specimen surface(s). In addition, 
power spectral density plots of input and response functions were recorded. 
This was conducted using a spectral dynamics Model 301c real time analyzer. 

A SD 302c ensemble averager was used in conjunction with the 301c to pro- 
vide a time-averaged value of power spectral density. The time-averaged 
value of power spectral density were used in evaluating the modal density 
of the specimen. In averaging, the discrete frequency components of the 
response function are enhanced while purely random components diminish. 
Resolution of data in the frequency domain was obtained by the use of a 
SD 107 low frequency translator. The translator provided a means of analyz- 
ing a narrow band of data at any frequency. 

After completion of tests to determine the response ratios and the 
modal density under narrow band random excitation, each specimen was sub- 
jected to sinusoidal input of constant force. The force level was main- 
tained by a SD 105 servo amplitude controller. A very stable SD 104-5 
oscillator was used as a function generator and a Hewlett Packard Model 
5323A electronic counter provided an accurate means of frequency readout. 

The counter displays the frequency to seven place accuracy in .4 sec. 

At the major modes of frequency response on the structure being excited, 
the 3 db bandwidth was obtained as a measure of damping. Several criteria 
were established for the selection of frequencies where damping was measured: 

1. The frequencies represented major modes of the surface being 
excited. 

2. No other modes existed close to the selected frequency. 

3. The 3 db down points were nearly equally spaced about the 
selected natural frequency. 
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A block diagram of the typical instrumentation is shown in Figure 7. 

For tests of the three simply connected structural models and the base 
plate, the specimens were supported at the four corners of the rectangular 
base plate. The supporting structure was non-resonant up to 2kc. However, 
some motion was detected at the attachment above this frequency. Between 
2kc and 6kc the amplitude of the attached corners was less than one tenth 
the input level. At frequencies above 6kc, corner motion reached unity 
with respect to the input excitation level. The overall effect of comer 
motion at high frequencies was not examined. However, it is expected that 
the high damping values measured above 6kc were in part attributed to the 
attachment. 

The vertical plate and the half cylindrical shell components were 
tested separately by suspending them from light strings. No problems 
related to the means of suspension were detected. 

In each test condition the shaker was in turn, attached through a 
.635 cm (.25 inch) diameter hole to each surface of the specimens. 

Models No. 1, No. 2, and No. 3 were first excited at the base plate with 
response ratios, spectral density measurements and damping ratios determined. 
The shaker was then suspended by shock cord and attached to the vertical 
plate, cylindrical half-shell, and open-end cylindrical shell, respectively, 
and a similar set of tests were repeated. 

Test Data of Model No. 1 

The first model (Figure 3a) has an irregularly shaped vertical plate 
welded to the base plate. The general configuration is similar to the 
model tested by Lyon and Eichler (Reference 4) except that the plate 
thicknesses are reduced. 

The measured values for loss factors (= 2c/c ) were averaged within 

r a 

each 1/3 octave frequency band. The apparent loss factor n ^ of the top 
plate, measured at the top plate of the two-plate system when the top plate 
is excited, is plotted in Figure 8 vs. 1/3 octave band center frequency. 

Also plotted in the figure is the dissipative loss factor n of the top 
plate component. Similarly, Figure 9 shows the plots of the apparent loss 
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FORCE LINK 



FIGURE 7 INSTRUMENTATION BLOCK DIAGRAM FOR MEASURING 

STEADY-STATE RESPONSE 
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factor n b of the base plate for the two-plate system and the dissipative 
loss factor of the base plate component. For most cases, the apparent 

loss factors are greater than the corresponding component loss factor. 

This is as it should be, as indicated by Equation (52) previously. The 
few data points which do not satisfy the above condition are considered 
unreliable and are rejected. For the steady-state response measurements 
the specimen was excited with 1/3 octave frequency band random force, 
and the rms response acceleration was recorded at selected locations on 
the specimen surfaces. Typical input PSD plot with center frequency at 
630 Hz is shown in Figure 10, while the corresponding response is shown 

in Figure 11. The ratios of the measured mean square values of average re- 

2 2 

sponse <V >/<V b > when the base plate was excited were plotted in Fig- 
ure 12 vs. center frequency of excitation (broken line). Also plotted in the 
figure are two sets of estimated response ratio data based on: 


1. Strong coupling assumption: using Equation (53) and measured 

, _ a 

loss factors n n t : 


< v t 2 > 


M b V - n b 

M 


(78) 


<Vb > "t t 

Weak coupling assumption: using equations (50), (56), and 

a 

measured loss factors n , n ^ • 

tb t 


<V > 

2 

<v/> 

b 


M, 


N 


M 


N, 


tb t 


tb 


(79) 


In Figure 12, a line corresponding to 


<v t > 

<v 2 > 

b 


M. 


M 



was plotted. This line serves as an upper bound of the actual <v c >/<v b > 
ratio for the weak coupling case in view of Equation (79). It may also be 
stated that the line represents a condition where the average modal energies 
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FIGURE 10 TYPICAL 1/3 OCTAVE 630 Hz CENTER FREQUENCY INPUT PSD AT LOCATION A 



Hz FILTER BANDWIDTH 



FIGURE 11 TYPICAL RESPONSE PSD AT POINT 2 DUE TO 630 Hz CENTER FREQUENCY INPUT AT A 


RESPONSE RATIO <V >/<V 



l 60 400 1000 2500 6300 16000 


THIRD OCTAVE BAND CENTER FREQUENCY Hz 

FIGURE 12 RESPONSE RATIOS OF THE TWO- PLATE SYSTEM 
WHEN BASE PLATE IS EXCITED 
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in the directly excited substructure and the indirectly excited substructure 
become equal. This condition makes it necessary that has to approach 

infinity in order for SEA to be applicable (see Equation (39) ) . Referring 
to Figure 12, the data seem to indicate that for the major portion of the 
frequency range, the computed response ratios based on the measured damp- 
ing factors of a weak coupling case are more close to the actual measured 
ratio, while the computed ratios based on the strong coupling case are 
higher except in the low frequency end. The difference in the computed 

ratio data is believed due to the fact that the lower thickness of the 

2 2 

top plate yields a low n t which gives rise to a high < v t > / < ^ rat *° 

based on Equation (78). As will be shown later, the trend is reversed 

(i.e., strong coupling equation gives better fit), when the top plate is 

excited. Our conclusion is that in this case, the top plate is virtually 

2 

cantilevered, the measured mean value < V > is sensitive to the locations 
where the measurements are made. In the high frequency end (e.g., 16,000 
Hz), the large deviation in the measured and computed results indicates 
the uncertainty in damping measurements based on which data Equations (78), 
(79) are applied. 

When the top plate was excited, the estimated response ratio 
2 2 

<v b >/ <V > may be obtained by using Equation (53) with the measured 
loss factors and p . In addition, instead of using measured 

values of (n - p^) , it may be evaluated analytically through the use of 
Equation (63) and the absorption coefficient Y. For the test model con- 
sidered, the values of Y may be computed based on Equation (77). The results 
of Y/Y q are plotted vs. one third octave band center frequency in Figure 13. 
As shown in the figure, the straight line result is based on the Kirchhoff 
plate theory and the curved line is based on the Mindlin theory. Similar 

to Figure 12, Figure 14 shows the measured and computed values of the 

2 2 

response ratio < >/<V > when the top plate was excited. The dotted 

line represents the response ratios computed based on the theoretical 
absorption coefficient Y as given in Figure 13. The computed response 
ratios based on Equations (78), (79) are also plotted in Figurel4. In 
this case, Equation (78), based on the strong coupling assumption, yields 
data which are closer to the measured response ratio data. 
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FIGURE 14 RESPONSE RATIOS OF THE TWO-PLATE SYSTEM WHEN TOP PLATE IS EXCITED 
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In order to evaluate the average mode-to-mode power flow coefficient 
g using Equations (47) and (50), the modal density is an important physical 
parameter. The modal density of a plate is given as: 

n p (f) = V3 (1 -v 2 ) A p /C p h p (80) 

where A p is the area of the plate, h p the plate thickness, v the Poisson's 

ratio, and C p the longitudinal wave speed of the plate material. Equation 

(80) gives n (f) = .232 mode/Hz for the base plate and n (f) = 0.1605 
P c 

mode/Hz for the vertical plate. The analytical modal density data may be 
compared with the PSD peak count data in a high resolution PSD plots such 
as those generated by the real time analyzer. Typical PSD plot of the 
vertical plate component covering the frequency range (580-680 Hz) is 
shown in Figure 15. The mode count yields 15 modes in 100 Hz band at the 
center frequency of 630 Hz, which compares well with the modal density 
data quoted above. 

It has been shown previously (Inequality (61)) that in order for the 
SEA method to be applicable in a connected structure, it is required that 
the two-way average power flow coefficients between the mode sets be equal 
and positive. This condition is now examined based on the measured and 
computed response ratio data. The values of the average mode-to-mode power 

flow coefficients g may be obtained as follows: 

8 bt . — using Equation (59) with measured data of response ratio 
2 2 

< v t > /< v b > when the base plate was excited, and n t . 

2* g tb — using Equation (59) with measured data of response 
2 2 

ratio <V^ > /< v t > when top plate was excited, and n^. 

— f 

§bt us ^ n 8 Equations (60) and (62) with the measured loss 

factors n t , and or using Equation (59) with computed 
data of response ratio < V ^ > as shown in Figure 12. 

^ • 8 t b using Equations (60) and (62) with measured loss factors 
h t * an d or using Equation (59) with computed data of 

2 2 

> /< V > as shown in Figure 14. 
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FIGURE 15 TYPICAL MODAL MAKEUP PLOT OF THE VERTICAL PLATE COMPONENT 



The results of and are plotted in Figure 16 while g^ and g ^ 
are plotted in Figure 17. As described above, is computed based on 

Equations (60) and (62). Referring to Equation (62), and g^ t are negative 
if the numerator on the right side of (62) is positive and the denominator is 
negative. In other words, g^ t is negative (and thus unacceptable in the 
physical sense) if the following is true: 


n 

t 




b> <o 


(81) 


For the two-plate system, it was found that in the frequency range (630- 
1600 Hz), condition (81) is true. As a result, g^ t data are not shown in 
Figure 16 corresponding to this frequency range. 


Test Data of Model No. 2 

The second model has a half cylindrical shell welded lengthwise to the 
base plate along the longitudinal boundaries of the shell (Figure 3b). The 
average values of the apparent loss factor n of the half shell of model 
Number 2, as well as the dissipative loss factor of the half shell com- 
ponent, are plotted in Figure 18. (Note that Figure 18 also includes one 
apnarent loss factor plot for model No. 3, to be discussed later in the 
report.) The average values of the apparent loss factor n ^ of the base 
plate of model No. 2 are plotted in Figure 19. 

The response ratio plots for model No. 2 are given in Figures 20 and 

2 2 

21. Figure 20 shows the values of < >/< v ^ > when the base plate 

was excited. Also plotted in the figure are the computed data based on 

Equation (53) and the measured loss factor data. The response ratio 

2 2 . 

plots < > /< > when the top shell was excited are given in Figure 

21, where one curve corresponds to the measured data, and the other curve 
corresponds to the comouted data based on Equation (53). 

The analytical expressions of the average input power per unit width 
in the plate and the average transmitted power in the half cylindrical shell 
are defined in Appendix IV by Equations (lV-4c) and (IV-34) respectively. 
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FIGURE 17 AVERAGE POWER FLOW COEFFICIENTS OF THE TWO- PLATE SYSTEM 
BASED ON COMPUTED RESPONSE RATIO DATA 
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FIGURE 19 LOSS FACTORS MEASURED AT THE BASE PLATE OF MODELS NO. 2 AND 3 
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RESPONSE RATIO <V S/<V, 



FIGURE 20 


RESPONSE RATIOS OF MODEL NO. 


2 WHEN BASE PLATE IS EXCITED 
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RESPONSE RATIO <V >/<V 
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FIGURE 21 RESPONSE RATIOS OF MODEL NO. 2 WHEN TOP SHELL IS EXCITED 
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Work was carried out to calculate the ratio of the average transmitted power 
to the average input power. For this computation, use was made of Equations 
(III-9) to ( III-13) , ( III-20) , and (III-21), corresponding to various input 
bending wave frequencies and angles of incidence for Model No. 2. The re- 
sults indicate a rather erratic variation of the average power based on 
strain energy rate in the shell with small values for some angles of inci- 
dence interspersed by very large peaks. The behavior is attributed to 
resonance effects in the cylindrical shell, most likely to the interaction 
of flexural and extensional waves. Because of the excessive oscillation 
of the power ratio as a function of the incidence angle, our conclusion is 
that it is impractical to compute the energy transfer coefficient of this 
model using the present analytical formulation. 


The modal density of the half cylindrical shell is given in Reference 


35 as : 


n c <£) = n p (f) -f 


■ P 

J ' jrrr ( 


de 


„ V 1 - <f r /£) 


s in^ 0 


(82) 


where 


n (f) 
P 


= mo 


dal density of an equivalent plate, Equation (80) 


f 


1 

2 7T R 


if 


(83) 


radius of the cylindrical shell 


e 


u/2 


f >f 


-1 


\ sin vm r f<f r 


The ratio n /n based on Equation (82) is shown as the solid line in Figure 22. 
° P 

The modal density function has a singularity at frequency f^. As an alternate 
approach the total number of modes below a frequency f is 


M (f) - n (f) f -~- 

P 71 


f V 


1 - ( f r /f) 2 sin 4 e de 


(84) 
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FIGURE 22 AVERAGE MODAL DENSITY OF THE HALF CYLINDRICAL SHELL OF MODEL NO. 



The average modal density n^ (f Q ) over a frequency bandwidth of A f with 
center frequency f may be obtained as follows: 


n 

c 




M(f^ 4* \ Af) - M(f Q - \ Af) 
Af 


(85) 


A computer program was prepared to compute the values of n^( f Q ) /n^( f Q ) f° r 
each 1/3 octave frequency band from f^ = 250 Hz to f^ = 16K Hz. The result 
is plotted in Figure 22 by the dotted lines. For the half cylindrical shell 
of model No. 2, f is equal to 3120 Hz. Based on the values of n^(f), the 
upper limit lines for the response ratios similar to the limit line of 

Figure 12 are drawn in Figures 20 and 21. The limiting line corresponds to 

2 2 2 2 
<V^ >/<V^ > = ^b^c^c^b F*-S ure 20 and corresponds to <V^ > /< V c > 

M c^b^ M b N c * n F *8 ure 21 • The average mode-to-mode power flow coefficients 

for model No. 2 are plotted in Figures 23 and 24. Similar to the two-plate 

system (Figures 16 and 17), the values of g , g^ are computed based on 

the measured ratios and the dissipative loss factors using Equation (59). 

— i —i 

The values of g^ c and g^ are computed based on the measured apparent and 
dissipative loss factors using Equations (60) and (62). 


Test Data of Model No. 3 

The third model features an open-end circular cylindrical shell welded 
to the base plate (Figure 3c). The average values of the apparent loss 
factor n cb of the cylindrical shell of the model were measured and plotted 
in Figure 18. Since no individual cylindrical shell component was fabri- 
cated, the values of dissipative loss factor n of the half cylindrical 
component of model No. 2 (with the same thickness) were used in the compu- 

tation. The average values of the apparent loss factor n of the base 

be 

plate of the model are plotted in Figure 19. 

The response ratio plots of the model are given in Figures 25 and 26. 

2 2 

Figure 25 shows the measured and computed values of < V > /< V, > when 

2 2 C b 

the base plate was excited. The plots of < > /< V c > when the cylindri- 

cal shell was excited are illustrated in Figure 26. 
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FIGURE 25 RESPONSE RATIOS OF MODEL NO. 3 WHEN BASE PLATE IS EXCITED 
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RESPONSE RATIO < V, >/<V 
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FIGURE 26 RESPONSE RATIOS TO MODEL NO. 3 WHEN TOP SHELL IS EXCITED 
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For a finite cylindrical shell of length Jt , Reference 55 gives the 
total number of natural frequencies below a given frequency f. The formula is 


M(f) = 


n (f) 

P TT 


% (£ > £ 


r 7 (2 u - 1) 


[-J-+ sin 1 (2 y - 1)J + Vu-U 2 | 


for f < f 


for f > f 


V = f/f 


( 86 ) 


where f is the reference frequency as defined by Equation (83), and n (f) 
is the modal density of a flat plate with the area (half the cylindrical 

surface). The average modal density n^ (f Q ) over a frequency bandwidth Af 
with center frequency of f^ may be obtained using Equation (85). With the 
modal density data, the upper limits of the response ratios for the sub- 
structures of model No. 3 are determined. These limits are plotted in 
Figures 25, 26. As shown in Figure 25, both the measured and computed 
response ratio data surpassed the upper limit curve for frequencies under 
2000 Hz. This fact can be interpreted by the observation that the inter- 
face between the plate and the cylindrical shell forms a boundary which has 
a dominating effect on the low frequency motions of the based plate. As 
a result, the SEA formulation in its present form is not applicable to the 
structural configuration in the low frequency region. 
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CONCLUSIONS 


In this report, the basic foundation and assumptions of the SEA method 
were studied in detail. A new formulation based on the strong coupling case 
has been generated and presented. As far as the power flow coefficient was 
concerned, the previous formulation based on the weak coupling case could 
be considered as a special case of the strong coupling case. On the other 
hand, in determining the kinetic energy transfer between connected sub- 
structures, the presentation of the coupling loss factors and the defini- 
tion of the apparent loss factors of a connected substructure were differ- 
ent between the weak and strong coupling cases. This deviation in turn 
caused a difference in the predicted substructure response ratios based 
on the experimental loss factors. 

In an effort to determine the power flow coefficients between two 
typical connected substructures, three simple models were made which 
included a two plate system, a plate and a half cylindrical shell system, 
and a plate and an end-connected cylindrical shell system. For the two 
plate system, the analysis was carried out to completion. The response 
ratio prediction based on the analytical formulation indicated that through 
the frequency range of interest, the analytical method was correct to 
within the order of magnitude of the mean square response ratios (see 
for instance, Figure 14). In evaluating the deviation between the 
experimental and analytical results, our conclusion is that part of 
the deviation was due to the limitation and arbitrary assumptions made 
in the Statistical Energy Analysis method. On the other hand, the measure- 
ment of the experimental responses was based on the averaging of the response 
data at a number of randomly chosen locations. Since the motion was not per- 
fected diffused, the experimental response data was dependent on the random 
locations selected. The existence of a partially diffused wave motion was con- 
sidered a contributing factor to the deviation between the experimental 
and analytical results. Our experience on the SEA method is that even 
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though the basic guidelines as developed in the report are satisfied, the 
predicted mean square values are usually only accurate to within the order 
of magnitude of the experimental results. 

For the plate and half cylindrical shell system, analytical formulation 
was carried out to define wave motions at various regions of the connected 
substructures. In the process, a set of definitive expressions was developed 
for the wave motions considering the interface. Based on the wave motion 
equations and the response formulation of Appendices III and IV, selected 
computation was carried out on the ratio of power input to the cylindrical 
shell at various frequencies and angles of incidence. Because of the exces- 
sive variation of the power ratio as a function of the incidence angle, it 
has been found impractical to compute the energy transfer coefficient using 
the present analytical formulation. On the other hand, for this model, the 
computed response ratios based on the substructure loss factors and the 
apparent loss factors compared favorably with the measured response ratio 
data (Figures 20, 21). As far as the average power flow coefficients were 
concerned, the variations between the experimental data and the computed 
data again covered a band whose maximum at any one frequency might differ 
from the minimum value by one order of magnitude. 

The third model consisted of a flat plate and a complete cylindrical 
shell with one end welded to the flat plate. For this structure, because 
of the strong interaction of the interface, it was found that the SEA 
method was not capable of reaching any intelligent prediction of the 
coupling loss factors, and the energy transfer coefficients. This fact 
was also confirmed by the guidelines which indicated that the SEA method 
in its present form was not applicable to the model under consideration. 

In exploring the application of the SEA method, certain limitations were 
encountered. In general, for a connected structure it is not uncommon that 
the predicted mean square response data and the power flow data differed from 
the corresponding experimental data by one order of magnitude. For this 
reason, the SEA method seemed at best to serve the purpose of trend predic- 
tion when the guidelines were carried out and presented in this report. 

These derivations illustrated the vibrational energy transfer mechanisms 
between the specific substructures under consideration. 
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APPENDIX I 


FLEXURAL WAVE PROPAGATION IN A TWO PLATE SYSTEM 
CONSIDERING ROTATORY INERTIA AND TRANSVERSE SHEAR 

Mind! In (Reference 50) deduced a two-dimensional theory of flexural 
motions of isotropic, elastic plate from the three-dimensional equation of 
elasticity. The theory includes the effects of rotatory inertia and trans- 
verse shear in the same manner as Timoshenko’s one-dimensional theory of 
bars. For this case, the three equations of motion may be expressed in 
terms of the plate-displacements as follows: 


D 

2 



-v) V 


+ ( 1 + v) 
T x 


d0 

3x 


* 3 Gh (V x 


+ d w \ ph 3 

dx/ 12 



5 I (1 - v) V% , +Cl+v) £0 

O T V . _ 


'- 1 '" (*y ♦ t;) 



* 2 Gh (7 3 w + 0) = ph 


at 2 


w 


where 


V = al 

dx 3 


9y 


(i-D 


( 1 - 2 ) 


(1-3) 


65 



H 

G 

D 

v 

h 

P 

w 

lb 1 b 

r x » r y 


a constant involving Poisson’s ratio (see Reference 50) 
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FIGURE 1-1 A TWO-PLATE SYSTEM 
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Let the incoming wave (see Figure 1-1) in Plate No, 1 be of the form 
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(I-7a) 


where, from equation (1-6), the wave speed c^. is given by 
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Associated with the normal displacements are the rotation angles and*l* ^ 


given by equations (1-1) and (1-3) as 
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Equation (I-7b) yields only one real wave speed, the value when the plus sign is 
used. In what follows, the range of to will be assumed to be governed by 
equation (1-9). Then the outgoing flexural waves in Plate 1 can be represented by 
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where from equations (1-1) to (1-3) and (1-6) 
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Note that the terms multiplied by C, are the so-called thickness-shear mode 
and are associated with a state of zero normal displacement. In the range of 
frequencies considered here only one outgoing wave is harmonic while the other 
two are of the Rayleigh type with an exponential decrease in amplitude as the 
distance from the joint increases. Similarly for Plate 2 the outgoing waves 
are given by 
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The nine constants of integration are determined in the present problem 
by the following conditions: 
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along x = z = 0. With the use of equations (1-5) and (1-16), equations (I-16d) 
and ( I— 16f ) may also be written as 
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The quantities w ] , and i|^ are the sum of the incoming and outgoing 

waves in Plate 1. The first of equations (I-16a) and equation (I-16b) state 
that Plate 2 is continuous so that lines along the y-axis originally perpendicular 
to the middle surface rotate the same amount. The first of equations (I-16c) 
and equation (I-16d) state that Plate 2 is supported along the y-axis in such a 
way that the normal deflection vanishes and the twisting moment is continuous. 

The remaining conditions imply that Plate 1 hypothetically extends to the 
middle surface of Plate 2 where it is fastened along the y-axis in such a way 
that a line originally normal to the middle surface has the same rotation in 
the x-z plane as a normal line in Plate 2 and remains in that plane. 
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Substitution of equations (I-7a), (1-8), ( I —10 ) , (I— 11), and (1-12) in equations 
(1-16) yields the following relations for the constants of integration 
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If the effects of the transverse shear deformation and rotatory inertia are 
omitted. Equation (I-17c) will be reduced to 


A 2 _B 2 = -2i(k 1 r/kg) cos 0 X [i C os 0 X - /l + sin 3 0 X 

+ r[i cos 02 - /l”+ sin 3 0 2 ]}’ L 
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Equation (I-17g), a special case of equation (I-17c), is identical to the 
result obtained by Lyon and Eichler in Reference 4. 

If the incoming wave is in, say, the z>0 side of Plate 2, the analysis 
is modi f ied by adding the incoming wave solution 
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to the outgoing wave solutions for the z >0 side of Plate 2 and deleting the 

incoming wave solution for Plate 1. Satisfaction of the boundary conditions 
then yields the following values of the integration constants: 
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If oa exceeds the value — = — in a plate, the wave speed given by 

h P (1) (2) 

Equation (l-7b) has two real solutions, say and . In addition, a 

(3) 

real wave speed can be associated with the thickness-shear mode. Then 

the outgoing waves can be represented by 
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Thus, different coefficient relations must be determined from the range 
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In the range of frequency value given by equation (I-22a) the waves in Plate 1 
have exponential decay and harmonic space variations while in Plate 2 the space 
variation is harmonic only* When the frequency is governed by Equation (I-22b) 
both plates have only waves with harmonic space variations* 
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APPENDIX II 


EFFECT OF FLEXURAL AND TANGENTIAL 
WAVE COUPLING ON THE TRANSFER FUNCTION 


In this study the effects of rotatory inertia and transverse deformation 
are neglected to better isolate the effects of flexural and tangential wave 
coupling. Consider a plate, the middle surface of which is the x-y plane. 

The pertinent equations of motion and force-strain relations are given by 

DV 4 w + ph -|^ = 0 (Il-la) 
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An incoming flexural wave can be expressed as 


w i = e 
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Outgoing flexural waves are given by 
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Outgoing tangential waves can be written as 
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FIGURE II-2 INTEGRALLY JOINED PERPENDICULAR PLATES 




Consider now an imcoming flexural wave in Plate 1 of the form 
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Outgoing wave displacements in Plate 1 are 
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For the x >o and x<o sides of Plate 2, the outgoing displacements are 
given by 
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with 
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At the joint the displacements and rotations must he 

con t inuous . Then 
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These eight deformation continuity conditions are supplemented by four force 
and moment continuity conditions (see Figure II-3) 



FIGURE II-3 JOINT FORCES 
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which can be expressed with the aid of the deformation continuity conditions as 
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The substitution of equations (II-5), (II-6), and (II— 7) into equations 
(II-9) and (11-11) yields twelve simultaneous equations for the determination 
of the integration constants. 
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The effect of flexural and tangential wave coupling depends on the 
parameter £ . The results for uncoupled motion are obtained with £ equal 
to zero. If £ is large, however, the effects of coupling can be signifi- 
cant. In order to get a better idea of the magnitude of £ , let it be 
rewritten as 
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so that £ will also be small if the mass ratio is of the order of unity. 
However, £ becomes large when is in the vicinity of unity. This is 
possible in the range of values of CO of interest for a certain range of 
values of For example, for a steel plate with 

E 2 = 30x10* psi = 20.69 x 10* N/cm 7 
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which is small again. 
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APPENDIX III 


TRANSMISSION OF A FLEXURAL WAVE IN AN INFINITE 
FLAT PLATE WITH AN INTEGRALLY ATTACHED HALF 
CIRCULAR CYLINDRICAL SHELL 

INTRODUCTION 

The structure under consideration consists of an infinite flat plate 
to which an infinite half circular cylindrical shell is integrally joined 
along the diametrically opposite generators (Figure III-l) . A harmonic 
flexural wave is produced in the plate and is scattered by the half cyl- 
inder, The object of the investigation is to determine the scattered 
wave distribution in the plate-cylinder combination. In what follows 
rotatory inertia and transverse shear deformations are neglected in both 
the plate and the cylinder. Flexural and extensional wave coupling is 
considered, however. The effects of transverse shear and rotatory inertia 
in the cylinder can be considered with the use of Reference 53, if neces- 
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FLAT PLATE ANALYSIS 


The plate is considered to be divided into three regions, region 1 
for which x >R, region 2 for which x<-R, and region 3 for which -R<x<R. 
In region 1 the incident and reflected wave deformations are given by the 
following equations (see Appendix II) : 
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In region 2, the outgoing wave deformations are given by 
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Finally in region 3, waves are directed into the interior from both 
boundaries so that the deflections may be written as 


89 



u° = P e “Ao;(x + R) + c o e Xco(x-R) + go e"^ ( * + R) 


( III-3a) 


v° = i 
P 


sin 0 


Zl I B° e “ Aw(x + R) - c° 


sm 0 


/c 


-o e ~Acj(x + R) _ e Mx-R)J e i«(t- * S ^ 0 ) 


( III-3b) 


A 0 e 


-ia>£2££ (x + R) 


+ D° e 


. cos 0 , 

0 (x-R) 


a 


A 0 e -CU C < x + R ) 


+ D° (x " R) 


^(t - 


( III-3c) 


with 


c 


w 



h 

P 


a 


yr 


sin 2 0 


A 


A 



(1-m 2 ) 

E_ 

E 

P 


P 


£ 


2(1 + j/ ) p 

EL — _E 

E 

P 


( III-4a) 

( III-4b) 
( III-4c) 

( III-4d) 


90 


CYLINDRICAL SHELL ANALYSIS 


The equations of motion of a cylindrical shell with transverse shear 
deformation and rotatory inertia neglected may be written as (Reference 54) 
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The forces needed are those on the straight generators of the cylindrical 
shell (Figure III-2) which are given by Reference 54. 




FIGURE III-2 STRESS - RESULTANTS ON GENERATORS 
OF CYLINDRICAL SHELL 
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Let the solution of Equations (III-5) be of the form 
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with U, V, and W constants. Then the substitution of Equations (III-8) 
into Equations (III-5) yields the following set of equations to be 
satisfied. 
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Thus, for Equations (III-9) to be satisfied y must be a solution of the 
biquartic equation 
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Equation (III-10) has eight solutions which are denoted by 
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Then the expressions for the displacements, slope change, and pertinent 
stress-resultants may be expressed as 
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(A e —Be I 
\ n n / 


y n 0 -Y n 0 

e + B e 
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* \ 3e n - 5 n 


The common factor e 


(* - H* 2 *) 


has been omitted from each of 


Equations (III-14) for simplicity. 
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BOUNDARY CONDITIONS 


The 24 constants of Integration in the preceding equations are deter- 
mined by the boundary conditions that must be satisfied at the two lines of 
contact of the plate and the half cylinder. These are (Figures III-l and 
III-3) 


U + : 

P 

= u° = 

P 

W 

c 


( III-15a,b) 

+ 

V : 

P 

= v° = 

P 
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c 


(III-15c,d) 
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W : 

P 

= w° = 

p 

V 

c 


(III-15e,f) 

+ 
9 x 

aw° 

= _R 

9x 

I 

R 

(%-■•) 

(III-15g,h) 

< 

+ 

1 

N a - 

0 

( III-15i) 

N + 

X 

- N° - 

X 

Qe = 

0 

( HI-15 j) 

N + 

xy 

- N° + 
xy 

\ = 

0 

( III-15k) 

M + 

X 

- M° + 

X 

M 0 = 

0 

(III-15/) 

9 — - 

- 7T/2 and 




U° 

P 
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P 

“ W 

c 


( III-16a, b) 

v° 

p 

= V = 

p 

u 

c 


( III-16c, d) 

w° 

p 

= w = 

p 

— v 

c 


(III-16e,f) 

aw° 

3x 

aw' 

= E. = 

a x 

■-*( 

a w c \ 

a e v c / 

(III-l6g,h) 
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0 


( I IT —161) 


Qx " <C + N 0 = 

N x ” N x " = 0 (III-16J) 

N° - n” - N fl . = 0 (III-16k) 

xy xy 0£ 

M° — M" - M fl = 0 (111-162) 

x x 6 

at x = -R, G = tt/ 2. From the force-strain relations (see Appendix II) 
and the continuity conditions for displacements, the last four of 
Equations (III-15) and of Equations (III-16) may be expressed as 


D p sf? (”p - ”p) - * 0 <m-i7a) 

-fc (■;-«;)- - o (m - l7W 

-T&nrj £- ft - ) + V « (III - 17c) 

D p al 5 ' ("p _ w p ) _ M e * 0 <m-l 7 d) 

at x = R, Q = — ^ and 

D THT (*° ~ w" ) - N fl =0 ( IH-18a) 

p ax 3 \ p p/ e 


af («;-«;)- 5 . - o (III - 18b) 

E 2 h 2 a / o - \ 

Wl +\) ~ ( P ~ V P / - = 0 

d p 5? ( b ; ■ "p ) + ° 


(III-18c) 
( III-18d) 
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<4 

(b) x = R, Q = - 

FIGURE III-3 JOINT FORCES 
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Substitution of Equations (III-l - III-3) and (III-14) into Equations 
(III-15a-h), (III-16a-h), (III-17) and (III-18) and manipulation of the 
resulting expressions then yields the following set of equations 


B + ± B _ = (< 


B' ± B = le“ 2 * wR ± y ) ( B ° ± C ° ) ± T^6 ( B ° ± C °) 


(III-19a,b) 


d ± B* 
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( III-19c, d) 
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A 
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. ^ d 



( A n " B n) cosh ( T n 


( III— 20a, b) 
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A 0 + 


D‘ = 1 L 


A n = 1 
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ico R cos 0 ^ _ e -2iwR cos0/cj f 


.( 1 + e -2i„ R c„sVc) (l + g l ( An + Bn)sinh(Ynf) 


K r 


2D 

p\ c 


4 

L 

n = 1 
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k + 


WRQ! 


[- 


y 2 e + v 4> 2 6 
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with 


A, 

A, 


- ( sln */c)’ 



(III-22a) 

X X 





i cos 0 



(III-22b) 
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L * 


I - 2 XcuR 

/ -2XcuR\ 


r sin0 / c 

\ l 6 

\ L + e / 




siny c 

(l~ + e - 2 ^ R )(i ±e - 2lwR ) 

(III-22c,d) 



X 



c 

a 

|l ^ e _ 2wotR/c 

j ± e '2i u>R cos <t>/c | 




— i cos 0 (l 

± e -2 w «R/c ) (l T e ~2ico R cos«/c) 

( III-22e, f ) 


It will be seen that Equations (III-20a, d, f, and h) constitute four equations 
for the coefficients (A r - B ) (n = 1, 2, 3, 4) while Equations (III-21b, d, 
f, and h) constitute four equations for the coefficients (A^ + B ) (n = 1, 2, 

3, 4). Thus, the solution of the problem has been reduced to the solution of 
two independent sets of four equations. 


105 



APPENDIX IV 


ENERGY TRANSFER IN VARIOUS STRUCTURAL SYSTEMS 


TIME AVERAGES OF HARMONIC COMPLEX QUANTITIES 


In what follows we shall be dealing with time averages of quantities 
of the form (Re A e ) (^e B where A and B are complex numbers. 
The integration of the quantities are carried out for one period T = — . 


/ 

icot \ 

/ icot 

\ 

1 

\ e 

Ae / 

\Ke Be 

) dt = 

4T 


/ „ tot 

~ - icjt N 

\ 



^ Be 

+ Be 

) dt 



+ ( A 

b+ab)I 

dt 

= 


4T * / 0 


Zt f 0 \ + 


A Be 


-2iwt 


AB + AB 


(IV-1) 


where the tilde over a quantity indicates the complex conjugate of that 
quantity. If A is equal to B, the right-hand-side of Equation (IV-1) 
becomes \ A A. 


AVERAGE INPUT POWER 


The average input power per unit width for a flat plate governed by a 
Kirchhoff plate theory is given by (Figure IV-19 and Appendix II) 



dt 



+ 


(2 - V) 


d 3 w 1 gw 


dt 


( IV— 2) 
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equation (IV-2) becomes 



( IV-4a) 


It should be noted that the input power is actually constant in this case. 
Since c is given by 


c 



( IV-4b) 


equation (IV-4a) may be written as 

P. = P ha) 2 c* ( IV-4c) 

1 


If shear deformations are included, the input power is given by (see 


Figure IV-lb and Appendix I) 



+ 


M 

xy 




) 


dt 


1 

T 



+ 



With the deformations given by (Appendix I) 


w = e 


-uo (. PcM icu(t+— ) 

~ y--—) e v c/ 


ib =0 


( IV-6a) 
( IV-6b) 
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Equation (IV-5) becomes 


■ PX (i- ¥) 


+ 

CJ 2 D 

o/d 

= 2^ 


Re e \ c / 

. ico/t + - ) 

Re le \ c / 

[>- *)*♦ sr] 


dt 
Phc' 


( IV-7) 


The wave speed c in this case is given by 


I = I 

c c* 


i + 




*)* 

+ h ( 
\ 

i + * 2 W 

} J 


z / 


( IV-8a) 


with 



$ = 


Then Equation (IV— 7) becomes 


( IV-8b) 


p h oPc* 


1 - 4* A ¥ «>1Vmi- ¥«vT 

yyx¥ 


^ (l - ^ (l + -^ * 2 )$ 


(IV-9) 

The quantity given is thus the ratio of the power required for the Mindlin 
and Kirchhoff plate theories to yield a flexural wave having a given trans- 
verse amplitude. The ratio is less than unity for all values of $, ranging 
from approximately 1 - | $ for every small $ to j — ■ for large It 

should be noted, however, that Equation (IV-9) is valid only for $ 

; I \ (Appendix I). 

•7 ( 1 - p )* 2 
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ORTHOGONALITY RELATIONSHIPS FOR AVERAGE POWER EXPRESSIONS 


We can show generally that the average power involved in the forces of 
one mode of motion moving through the displacements of another mode is zero. 
For the Kirchhoff flat plate, for example, with 


w = Re F(x) 
o 


ia ,( t _ Z_2±££) 

x) e \ c f 


(IV-10) 


with the average output flexural power given by an equation of the form of 
Equation (IV-2), and with use of Equation (IV-1) we have 


BK 


ia>D 
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ia?D r 
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Td 2 F 

y /co sin0\ 2 ~ 
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dx 
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«-»)(“ s ‘ n? ) 

dF 


dx 

[ d 3 F 

v / co sin0 \ 2 

dF 1 

L dx 3 

- c ) 

dx J 

f d 3 F 

~ d 3 F 

d 2 ] 

L ^ 

F " F - 

die 


d 2 f 

dx 2 


£ - (£ ? - i f )] 


(IV-ll) 


However, F and F satisfy the equation (Appendix II) 

(1 - sin 4 0) (F,F) = 0 ( IV-1 2) 


dV.F) 

dx 4 


- 2 


/ co sin0 \ 2 d ? (F.F) 
\ c / di? 


CO 
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If we multiply the equation for F by F, the equation for F by F, and sub- 
tract the two we have, after some manipulation 

d^F d 3 F _ d^F dF d 2 F dF 

dx L dx dx 3 dx 2 dx dx 2 dx 


dF p 


f * 


which states, then, that the average power is a constant and thus is inde- 
pendent of x. The implication of this result is that any product terms in 
Equation (IV-11) which are functions of x can be ignored. This further 
implies that the modes are orthogonal and that only exponential solutions 
having pure imaginary characteristics contribute to the average power. 

For the Hindlin plate the deformations are given by (Appendix I) 


w 

= Re W(x) 

icj 1 
e 

(t-: 

Y sin0 
c 

j ( IV— 14a) 

ip 

T x 

= Re X(x) 

ic o | 
e 

(t-: 

i sin0j 

( IV— 14b) 

ip 

y 

= Re Y(x) 

icol 
e 1 

I *- 1 

r sin0 j 

(IV- 14c) 


with the average output flexural energy given by an equation of the form 
of Equation (IV-5), the use of Equation (IV-1) yields 
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and X, Y, and W satisfy the complex conjugates of Equations (IV-16). Now 
multiply the three equations of the first set respectively by X, Y, and W, 
add the three resulting expressions, and subtract from it the sum of the 
three equations of the second set multiplied by X, Y, and W respectively 
to obtain the result that the derivative of Equation (IV-15) is zero. 

Thus, conclusions similar to those for the Kirchhoff plate hold. 


112 



For inplane motion of either a Kirchhoff or Mindlin plate the average 
power involved in extensional motion is given by (Appendix II) 

p = K /" / au dv \ du . il£ / du + av \ dv Hf . 

r E t J \dx d y ) at 2 y ay ax ) at 


with u and v given by (Appendix II) 


u- ReU(x) 

v = Re V(x) M- *■*%*) 


Equation (IV-17) becomes 


n iKcj ) / dU 

P E ' — (to + iV 


ivULSiSl V u - ( 4“ 

c / \ dx 


_ jy-feLsinp v 
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l-V \( dV 
2 \ dx 


+ i u v 


cj sm 




iKoj dU 

4 dx 


- »r 

dx 


1 + V 
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. co sin0 


+ UV 


¥(' 


- v — 

dx 


( IV-19) 


But U and V satisfy the equations (Appendix II) 


TT 1 + V . Cl) sin0 dV 
u — 1 d^ 


(IV-20a) 
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1 + V . w sin0 du , 1 - V d 2 V _ ,.,2 f sin 2 0 _ (1- ^)P | v = n 

' 2 1 “c dx + 2 dT 7 IfT 3 E J 

(IV— 20b) 

while U and V satisfy the complex conjugates of Equation (IV-20). 

Then multiplying Equations (IV-20) by U and V respectively, adding the 
two equations and subtracting the sum of the complex conjugates multiplied 
by U and V respectively yields the result that the derivative of Equation 
(IV-19) is zero. Thus, the same conclusions apply as for the Kirchhoff and 
Mindlin plate average bending power. 

Finally the case of the integrally connected half cylindrical shell 
and flat plate combination can be shown to yield similar results. The 
average power per unit length of the half cylinder is given by (Appendix III) 



( IV-21) 
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With u^ f anc j „ given by (Appendix III) 


u c = U(0) e 


; i (wt -4>0 

v c - V<9) e 1 ** 
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Equation (IV-21) becomes, with the use of Equation (IV-1), 
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( IV-22a) 
(IV-22b) 
( IV-22c) 
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(IV-23) 


Now U, V, and W satisfy the Equations (Appendix III) 
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(IV-24b) 
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| V c + k } U + 

(i+ 2 k#) 


dV 

d0 


+ k~7 + 2 (1 ->P 2 ) k + [l + k(l + 1 P 4 ) — S2 2 ] W = 0 

( IV-24c) 

and U, V and W satisfy the complex conjugates of Equations (IV-24). Multiply 
Equation (IV-24c) by W and subtract from it the sum of Equation (IV-24a) mul- 
tiplied by U and (IV-24b) multiplied by V. Subtract the complex conjugate 
of the resulting expression from it to obtain the result that the derivative 
of Equation (IV-23) with respect to 6 vanishes. The usual conclusions follow. 


AVERAGE TRANSMITTED POWER 


The results of the preceding section can be applied to obtain expressions 
for the average output power in the various structural systems under consideration. 

a) Perpendicular Kirchhoff Plates with no Motion of the Joint Line 

For this problem the outgoing wave displacement function is given by 
(Reference 49). 
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o 
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sin0 o \ 

cf I 


(IV-25a) 
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for x >0 and by the negative of the expression for x< 0, with 




a 2 = + sin 0 2 


( IV-25b) 


'* = 




( IV-25c) 


/ P, h 

2 2 


The average outgoing power is given by Equation (IV-2), for x>0 and by 
the negative of Equation (IV-2) for x <0. Then Equations (IV-11) and (IV-13) 

yield the total output power in both parts of plate 2 as 


2 P ? h 2 co 2 c* C 2 C 2 cos 0 2 


C >V / 

if sin0 2 = 2/c* sintf^cl 


(IV-26) 


if sin0 2 = C 2 /c* 


sin0 > 1 
l ? 


b) Perpendicular Kirchhoff Plates with Motion of the Joint Line 


For this problem the outgoing wave displacements are (Appendix II) 
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B - e V x 


sm^c* 2 / 2 J 
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[(a; cos VcJ + a, e a,a - x '' c ? ) e lu, ( t - )j 


( IV-28b) 


(IV-28c) 


The average output power is given by the sum of Equations (IV-2) and (IV-17) 
for x>0 and by the negative of these expressions for x< 0. Then the use of 
Equations (IV-11) and (IV-19) and the conclusions reached in the preceding 
section yield 

P o --(P B + P S ) (IV-29a) 

where 
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sin 2 0 


(IV- 29 c) 

c) Perpendicular Mindlin Plates with no Motion of the Joint Line 


The average outgoing power per unit plate width for a plate edge with 

x >0 is given by Equation (IV-5) and by the negative of Equation (IV-5) for 

a plate edge with x<0. The deformations are given by Appendix I, in the 

C 2 

case when sin0 2 = — sin0 < 1, as 
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(IV-30c> 




(IV-31a) 





( IV-3 lc) 



Then it can be shown from Equations (IV— 15) and the conclusions of the pre 
ceding section that 

P q = -2 P 2 h 2 (J 2 c* cos0 2 A + 2 A+ • 

1 -^( 1 - 





* / A A - B B \ (IV— 34) 

\ n n n n / 


where the summation is over all values of n for which y ' is pure imaginary. 
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